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Abstract 



> 

, We present a reciprocal-space pseudopotential scheme for calculating X-ray absorption near-edge 

• structure (XANES) spectra. The scheme incorporates a recursive method to compute absorption cross 

| section as a continued fraction. The continued fraction formulation of absorption is advantageous in 

£ — , that it permits the treatment of core-hole interaction through large supercells (hundreds of atoms). 

The method is compared with recently developed Bethe-Salpeter approach. The method is applied 
CN| , to the carbon K-edge in diamond and to the silicon and oxygen if-edges in a-quartz for which 

polarized XANES spectra were measured. Core-hole effects are investigated by varying the size of 
the supercell, thus leading to information similar to that obtained from cluster size analysis usually 
performed within multiple scattering calculations. 



1 Introduction 

Ch . 

XANES spectroscopy is a powerful technique to probe the empty states in solids. Moreover XANES 
gives information about the atomic arrangement around the probed atom up to the medium range 
order (i.e., about 8 A). Nevertheless the interpretation of XANES spectra is not straightforward and 
often requires sophisticated simulation tools. To calculate a XANES spectrum, one needs to compute 
rS \ the absorption cross section. This is given by the Fermi golden rule as a sum of probabilities per unit 

of time of making a transition from an initial state to an unoccupied final state through an interaction 
Hamiltonian [Q]. Electronic transitions involved in X-ray absorption spectroscopy are mainly governed 
by the electric dipole operator. The main difficulty of any absorption cross section calculation lies in 
the solution of the Schrodinger equation for the final and initial states. The choice of the method used 
to solve the Schrodinger equation depends on the localized or delocalized character of final states. 
Crystal field multiplet theory || is adapted to the calculation of localized final states, i.e., with 
strong electron-electron interactions. It is typically the case of Z/2,3 edges of transition elements || 
and M.4,5 edges of rare earth elements If the final state is delocalized (the case of K and L\ edges) 
multielectronic interactions are weak and monoelectronic approaches based on the density functional 
theory (DFT) are usually employed. Among monoelectronic methods, one distinguishes the real 
space (cluster) approach from the reciprocal space (band structure) approach. Real space multiple 
scattering theory ^ |(| Q has been extensively used in the past twenty years. However, multiple 
scattering theory suffers from an unavoidable drawback, i.e., the "muffin-tin" approximation of the 
potential. Recently, the "muffin-tin" approximation has been overtaken in a real space approach by 
using the finite difference method (MDF) to solve the Schrodinger equation Although it yields 
good results B, the method presented in Ref. [^] requires significant computing power. At present 
this is what restricts the applications to small clusters (around 50 atoms without any symmetry, at 
maximum) . 

With regards to band structure calculations, local projected densities of empty states (commonly 
called LDOS) have been often used to interpret XANES spectra (see for instance Refs. Jl^, [y], |l^]). 
Recently several band structure codes have incorporated the radial matrix elements [Il3| to properly 
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calculate XANES or ELNES (electron energy-loss near-edge structure). It is notably the case of 
Wien2k jlj], which is based on the full-potential linearized augmented plane wave method the 
orthogonalized-linear-combination-of-atomic-orbital program developed by Mo and Ching ]16| , [nj 
and CASTEP, whose ELNES part has been developed by Pickard 63, [H], po| |. This last approach 
uses pseudopotentials and reconstructs all electron wave functions within the projector augmented 
wave (PAW) method of Blochl Within these ab initio schemes the treatment of core hole- 

electron interactions is usually achieved through the inclusion of an excited atom into a supercell. 
The absorbing atom is then considered as an impurity. Beside these methods, Shirley j22j and 
more recently Soininen and Shirley |%| have proposed alternative formalisms that treat core-hole 
interactions in a different way. Indeed their methods incorporate core-hole interactions in the two- 
particle Bethe-Salpeter (BS) equation. Such methods are of great importance for optical absorption 
calculations ^5|, ^6|, However, in the case of K-edge XANES spectra, the core- hole is frozen 
at one atomic site and the two-particle BS approach can be reduced to a single particle calculation. 

By using supercell approaches, large systems (several hundreds of atoms) can be a priori treated. 
The charge density for such large cells is easy to calculate self-consistently by using pseudopotentials. 
On the other hand, the computation of the cross section is at present limited by the diagonalisation of 
the Hamiltonian for several empty states at many fc-points in the Brillouin zone. Indeed in Refs.Jl^, 
[lil pseudopotential calculations of ELNES spectra in diamond and c-BN have been performed with a 
small 16-atom supercell. Here we propose to use the recursion method of Haydock, Heine and Kelly 
|jo| , |3l| , to make the calculation of the cross section computationally tractable for larger supercells. 
Thanks to this iterative technique, the calculation time of XANES becomes negligible compared with 
the calculation of the self-consistent charge density in large systems. 

We incorporated the recursion method into a pseudopotential scheme. We first show in section 2 
that the absorption cross section can be easily computed within a pseudopotential formalism using 
the PAW method. Then, by using the recursion method, we express the absorption cross section 
as a continued fraction. We present applications to K-edges in section 3. Two materials have been 
chosen : diamond which is an interesting example for method testing, and a-quartz for which we 
have measured polarized XANES spectra. 

2 Method 

2.1 X-ray absorption cross section in the impurity model 

In a monoelectronic approach, the X-ray absorption cross section o~(uj) can be written as 



o-(uj) = 4n 2 aohuY^ \Mi^f\ 2 5(Ef - E t - Huj), (1) 
/ 

where hu) is the photon energy, ao is the fine structure constant and Mi-,; are the transition ampli- 
tudes between an initial core state with energy Ei, localized on the absorbing atom site Ro, and 
an all electron final state \ipf) with energy Ef 



Afi_/ = {i>f\o\i>i). (2) 

D is a transition operator coupling initial and final states. In the electric quadrupole approximation, 
D is given by £ ■ r(l + |fe ■ r), where i and k are the polarisation vector and the wave vector of the 
photon beam, respectively. Within the frozen core approximation, is a core state that can be 
taken from an all electron ground state atomic calculation. In the impurity model, \tpf) is an excited 
empty state that is solution of the Schrodinger equation for a potential that includes a core-hole on 
the absorbing atom. 



2.2 PAW formalism 

In the following we want to ex pre ss the transition amplitude (EqJ^) within the PAW formalism, as 
originally described by Blochl [|2l). We only recall the main aspects of the method that are needed 
to give a simple expression for the Mi_>/ terms. 

Within the PAW formalism, the final state all electron wave functions \ipf) are related to the 
corresponding final pseudo wave functions \ipf) through a linear operator T: 

m=r\i> f ). (3) 
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T differs from identity by a sum of local atom-centered contributions, that act only within spherical 
core regions centered on each atomic site R, called augmentation regions or SIr: 

T=l+Y, [|0K,»> - |£*,n)] (PH,»|. (4) 
R,n 

Here \4>R,n) and \4>R,n) are the all electron and pseudo partial waves, respectively, which coincide 
outside Qr. The vectors {pr,„\, called projector functions [plj, are equal to zero outside fi^ and 
satisfy the condition (j>n,n\4>R' ,n') = SRR'Snn' ■ The index n refers to the angular momentum quantum 
numbers (£, m) and to an additional number, used if there is more than one projector per angular 
momentum channel. The \4>R,n) form a complete basis for any physical non-core all electron wave 
function within Therefore the \tf>R,n) are also a complete basis for any physical pseudo wave 

function \ip) within Qr, i.e. for any (r\xR) function centered on an atomic site R and equal to zero 
outside Qr, 

J2WpR,n)((f>R,n\XR) = WXR)- (5) 

n 

Substituting Eq. (^) in Eq. Q and then Eq. (||) in Eq. (0), the transition amplitude 
becomes 

Mi->f = $ f \Q\ipi) 

+ ^$f\pR,n){(t>R, n \0\'>Pi) 
R.n 

-^2$f\pR,n)$R,n\&\i ) i)- (6) 

R,n 

In Eq.| the initial wave function (r\ipi) is localized on the site of the absorbing atom, R , then only 
the Ro term has to be considered in each sum. Furthemore it should be noticed that (r" |£5 is zero 
outside the O.r region. Therefore we can make use of Eq.(^) for the third term of Eq.(g), which 
thus vanishes with the first term. The transition amplitude is then reduced to one term. Now, 

introducing 

\<PR ) = Yl \PR ,n){<f>R ,n\£>\i>i), (7) 

n 

we obtain the following simple expression for the X-ray absorption cross section 



<t(u) = 4TT 2 a huJ2 Mf\&Ro)\ 2 5{Ef -Ei- huo). (8) 
/ 

The calculation of XANES spectra from Eq.^ presents us with the problem of determining many 
empty states |V>/), as has been mentioned in section 1. Indeed, the addition of unoccupied bands 
significantly increases computing time and then limits the size of the supercells. In the following, 
we show how the recursion method permits the cross section (Eq.^) to be rewritten as a continued 
fraction, so that only occupied bands have to be calculated. 

2.3 Recursion method 

The recursion method comprises of a powerful recursive algorithm that can be applied to a Hermitian 
matrix in order to transform it into a tridiagonal form. This method has been widely used in solid 
states physics. The pioneering works in this field were mainly by Haydock, Heine and Kelly [u). 
The first application to absorption cross section calculation was carried out by Filipponi within the 
multiple scattering formalism ^3]. In a close field Benedict and Shirley have implemented it for 
E2(<^) calculation including core-hole interaction j3^] within the BS approach. For simplicity, in the 
following, we assume that the norm of all electron partial waves coincide with that of the pseudo- 
partial waves. 

In order to use the recursion method in our scheme, we shall introduce the Green operator into 
Eq.(|^) by making the substitution 

$f)S(Ef -Ei- fuu)(4> f \ = --3m[G(E)], (9) 
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with 



G(E) = {E-H + i-y)~\ (10) 

In Eq.© G(E) is the Green operator associated with the pseudo-Hamiltonian H — HT , which 
is Hermitian. The energy E is given by E = Ei + hu> and 7 is an infinitesimal positive number. The 
cross section (Eq.^) can be rewritten as 

g(uj) = -47rao7iw3m [(£r„|(.E - H + i7) _1 |£R }l • (H) 

Following the original work of Lanczos J34|, |3^], the recursion method sets up a new basis in 
which the pseudo-Hamiltonian H, has a tridiagonal representation, from which the matrix ele- 
ments {ifR \(H — E — i"/)~ 1 \<fR ) are very simply derived. This new basis (called Lanczos ba- 
sis in the following) is obtained by the repeated action of H onto the normalized initial vector 
|uo) = |^h )/\/ (<?H l<?-Ro) through the symmetric three-term recurrence relation 

H\m) = a,i\ui) + 6<+i|«i+i) + 6i|ui_i). (12) 

where {ai} and {bi} are two sets of real parameters given by dj = (ui\H\ui) and bi = (iij|H|ltj-i) = 
(ui-i\H\Ui). The tridiagonal matrix representation of H in the {|ltj)} basis then leads to the following 
form of the matrix elements of Eq. ([ll ) 

@ Ro \(H -E- ijy^ Ro ) = ( ^ J ^ o) b2 • (13) 

a - E - i-y 1 j-s- 



A simple terminator has been used to finish the continued fraction |33j]. In particular, if TV is the 
number of iterations required to converge the calculation, we consider that the coefficients (a,i,bi) 
are equal to (a^f,&jv) for i > N, this leads to an analytical form of the terminator. It should be 
noted that the number of iterations iV strongly depends on the broadening parameter 7. With the 
iterative technique of Haydock the main part of XANES calculation involves the computation of the 
Hamiltonian acting on a single vector. This means that the computing time is considerably reduced 
compared to that can be required by Eq. ^ with an explicit diagonalisation. 

2.4 Comparison with Bethe-Salpeter approach 

Recently, the BS approach has been successfully applied to the calculation of XANES |^2|, |23| . Since 
the formalism of such two particles excited states calculations has been recently detailed in a review 
1^1, here we only recall the main aspects. The BS procedure consists of a ground state DFT 
calculation in order to obtain the Kohn-Sham orbitals, a GW calculation to correct the eigenvalues 
|^?| , and a solution of BS equation using GW eigenvalues and a static ally screened potential. 

Our method differs from the BS approach of Soininen and Shirley p3] by three points. The first 
two differences are substantial, the last one is just methodological. In particular, (i) we use the DFT 
eigenvalues, whereas Ref. |23| uses GW corrected eigenvalues JS?]]. If the GW correction on the 
empty states can be described by a rigid shift, as it is often the case, this difference will not affect the 
XANES spectra, (ii) We screen the core-hole by the valence electron response at all orders computed 
self-consistently within DFT. In Ref. |^|, the core-hole is screened at the linear order using the 
random phase approximation |^| . (iii) We develop the final state wave functions \ipf) in a plane 
wave basis set of a supercell. In Ref. |p3| , the BS calculation uses as basis the eigenstates of the 
ground state Hamiltonian. 

In the next section, our results will be compared with BS calculations as far as possible. 



3 Applications 

Our method was applied to the carbon K-edge in the diamond phase and to the silicon and oxygen 
if -edges in a-quartz. XANES calculations were performed in two steps: (i) the evaluation of the 
self-consistent charge density for supercells including one Is core-hole, (ii) the construction of a 
converged {|ui)} basis leading to the tridiagonal representation of the pseudo-Hamiltonian. The 
transition matrix elements were determined in the electric dipole approximation. Calculations were 
carried out using an ab initio total energy code based on density functional theory within the local 
density approximation (LDA)|^]. Norm-conserving Troullier-Martins pseudopotentials |^] with a 
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single component for each (I, m) component were used. For the excited atom a pseudopotential with 
only one Is electron was generated. In both diamond and a-quartz cases, the size of the supercell 
was increased until the neighboring excited potentials did not interact with each other. The wave 
functions were expanded in plane waves with an energy cutoff of 50 Ry and 70 Ry for diamond and 
a-quartz respectively. Reciprocal-space integrations were performed using a Monkhorst-Pack fc-point 
grid [0. 

3.1 Diamond 

Figure |l| shows calculations performed with four different sizes of supercell, all including a Is core-hole 
on one of the carbon atoms. The size of the primitive cell - trigonal cell (a=2.524 A) containing two 
carbon atoms - was successively multiplied in the three directions by a factor two (2x2x2 supercell, 
16 atoms), three (3x3x3 supercell, 54 atoms), four (4x4x4 supercell, 128 atoms) and five (5x5x5 
supercell, 250 atoms). Further increasing the supercell size does not provide significant modifications 
in the calculated spectrum, fc-point convergence of the charge density was obtained with a 2x2x2 
fe-point grid for both 2x2x2 and 3x3x3 supercells, and with only one fe-point for larger supercells. 
fe-point convergence of the Lanczos basis was reached with a 10x10x10 fc-point grid for the 2x2x2 
supercell, a 8x8x8 fc-point grid for the 3x3x3 supercell, and 6x6x6 fc-point grids for both 4x4x4 
and 5x5x5 supercells. The number of iterations of the Lanczos basis was checked with 7 =0.3 eV. 
The required number of \ui) vectors varied from around 800 to 1400, while increasing the supercell 
size from 16 atoms to 250 atoms. 

Calculations are compared with the carbon AT-XANES spectrum of diamond measured by Ma 
et al. j42|. The overall agreement between the experimental spectrum and the 5x5x5 supercell 
calculated one is quite good. In particular, the features from d to i are correctly reproduced as 
well in relative intensity as in energy position. Nevertheless in the first four eV of the spectrum the 
agreement between experiment and calculation remains unsatisfactory. The intensity of peak b is 
exagerated and the exciton peak a is not reproduced. Looking at the various supercell calculations, 
one notices that the first 25 eV of the spectrum are affected by the size of the cell, i.e. the interaction 
between core-holes belonging to neighboring cells. This gives information about the volume seen 
by the photoelectron during the absorption process: the photoelectron interacts with atoms located 
within a 6 A radius sphere centered on the absorbing atom. 

Several calculations of the carbon A"-edge in diamond have been reported in the literature. Weng 
et al. jl3| present pseudo-atomic-orbital density of p empty states together with multiple scattering 
calculations. On the one hand, "muffin-tin" potentials as required in multiple scattering theory 
appear to be inappropriate in that case. This result was confirmed by a 300 atom cluster calculation 
performed by the FEFF8 code |?], E3. On the other hand, although it neglects core-hole effects, 
LDOS calculation by Weng et al. [ |43[ seems to correctly reproduce features h and i. Pickard]!^] 
shows that including core-hole effects into a 2x2x2 supercell significantly improves the agreement 
with experiment. However, as the supercell size increases, peak b becomes much too intense (figure 
|l|). This may suggest an over estimation of the core-hole effects in that case. This result has also been 
observed by Soininen and Shirley^] while core-hole interactions is treated within BS approach. Their 
carbon K edge calculation in diamond, is however very similar to our 128 and 250 atom supercell 
spectra. It seems that the corrections brought to excited states by GW calculations has only a 
negligible influence on the energy positions of calculated features. This result is not surprising since 
previous studies^, ^(| have shown that the GW approximation used for band structure calculation 
in diamond, amounts to applying a scissors-operator shift which is valid to within a 0.5 eV error. This 
error is comparable to the energy resolution of XANES features. The exageration of peak b in both 
calculations requires more understanding and the problem of the treatment of core-hole interaction 
in diamond is at this time still unclear. 



3.2 a-quartz 

In the case of a-quartz, calculations are compared with polarized Si and O A"-XANES experiments, 
a-quartz (a-Si02) has a trigonal space group with a hexagonal unit cell. Since the point group 
is 32 (or D3), a-quartz is a dichroic compound Q. Therefore any absorption spectrum can be 
expressed as a linear combination of the two spectra <xii and a± corresponding to the polarisation 
vector e of the photon beam parallel and perpendicular to the ternary axis of the hexagonal cell (i.e. 
the c axis), respectively. While calculations of isotropic XANES spectra of a-quartz have been the 
subject of several studies in the past [ |47| , p^, |i"9] , pc| |Tij ] , little attention has been paid to the angular 
dependence of absorption. Lagarde et al. plfj show that the differences observed between crn and a± 
spectra are similar to that observed between the spectra of silica and densified silica. In Ref. [5lL p2] 
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polarization effects are interpreted in terms of geometrical considerations. However, to our knowledge 
no convincing calculations of polarized data have been reported yet. 

Si and O K-edge measurements were carried out respectively on the SA32 and SA72 beam lines 
at the Super- ACO facility of LURE (Orsay, France). The storage ring was operating at 800 MeV 
(A c = 18.6 A) and 200 mA electron current. The sample was a synthetic (1010) oriented single crystal 
of a-quartz. The SA32 beam line was equipped with a double InSb (111) crystal monochromator. 
The instrumental resolution was estimated to be around 0.7 eV. The incident beam was monitored by 
an ionisation chamber while the absorption of the sample was measured in the electron drain current 
mode. Energy calibration was checked after each scan by measuring the c-Si spectrum for which the 
inflection point of the absorption edge was set to 1839 eV. The spectra were collected over a photon 
energy range of 1830-1910 eV with 0.2 eV steps and 1 s integration time. The sample was placed on 
a rotating holder, the rotation axis of which was parallel to the photon beam. The (lOlO) plane of 
the sample was set normal to the photon beam and the spectra were measured for the polarisation 
vector e either parallel to the c axis (<tm) or perpendicular to the c axis We verified that the 

combination |ctx + reproduced the powder sample well. The X-ray natural linear dichroism 
(XNLD) is obtained by the difference cxii — a±. The SA72 beam line was equipped with a toroidal 
grating monochromator. The estimated resolution was around 0.3 eV. Polarized O Tf-XANES were 
measured in the fluorescence mode over a photon energy range of 525-580 eV with 0.11 eV steps. 
The first 12 eV were slightly affected by self-absorption and the signal to noise was not good enough 
to permit an unambiguous extraction of the XNLD. 

Calculations were performed with various sizes of cell built from the cristallographic parameters 
given in Ref. ^3|. Convergence in terms of supercell size was obtained for a 2x2x2 supercell (72 
atoms) at both silicon and oxygen 7f -edges. In order to show the influence of core-hole effects we also 
present results obtained from the primitive cell (9 atoms) with or without core-hole (ground state 
calculation). Charge density was calculated with 2x2x2 fc-point grid for primitive cells (with or 
without core-hole) and only one fc-point for the 2x2x2 supercell. The Lanczos basis was determined 
with 4x4x4 fc-point grid for the primitive cell and 3x3x3 fc-point grid for the supercell. The 
broadening parameter 7, was set to 1 eV, this led to the construction of around 500 and 600 \m) 
vectors for the primitive cell and the supercell calculations respectively. 

Experimental Si K-edge X-ray absorption spectra are compared with 2x2x2 supercell calculations 
in figure^. From top to bottom one can see cry contribution (figure ^a), a±_ contribution (figure^)), 
and the XNLD signal (figure |^c). Figure |^c shows that a-quartz presents a strong angular dependence 
at the silicon K-edge, close to 9% of the white line's intensity. Experimental and calculated spectra 
- polarized and dichroic- are in excellent agreement. The strength and the energy position of all 
the features are correctly reproduced. At the oxygen A'-edge (figures ^a and ^b) the agreement with 
experiment is also satisfactory, au and a± spectra differ essentially in the 12-19 eV energy range. 
These differences are well reproduced by the calculations. 

Figure ^ compares three calculated a± spectra at both silicon and oxygen isT-edges: the 2x2x2 
supercell spectrum (solid line), the primitive cell calculation including core-hole effects (dashed line) 
and the ground state calculation (dotted line). At the silicon K-edge (upper panel), the ground state 
calculation does not reproduce the intense white line, which is clearly a signature of core hole-electron 
interaction. Including a Is core-hole on one of the three silicon atoms of the unit cell drastically 
modifies the spectrum and results in the appearance of the white line. Increasing the size of the cell 
yields changes that are essentially in the 8-18 eV energy range, showing that the resonances occuring 
in that region are related to medium range organisation around the absorbing atom. At the oxygen 
K-edge (lower panel), although core-hole effects are less important than at the silicon K-edge, they 
definitely have to be taken into account in the calculation. Furthemore our simulations show the 
sensitivity to medium range order of the 12-20 eV region just above the main peak. 

The importance of core-hole effects in a-quartz was first noticed by Jollet and Noguera ^7| . This 
point was emphasized in multiple scattering studies Q . Our calculations confirm these previous 
considerations and demonstrate that reciprocal space calculations can also evidence the influence 
of medium range order, as is usually shown in cluster size analysis performed within the multiple 
scattering formalism M|. 

Figures |^ and ^ display a good agreement between calculated and experimental spectra. Nev- 
ertheless the agreement in energy position may be discussed to a certain extent. Indeed calculated 
features appear slightly too contracted compared with experiment. For example, 1 eV shift is ob- 
served at 25 eV in figure ^. This contraction becomes more pronounced as energy increases. This 
small discrepancy could be due to the use of DFT eigenvalues. Unfortunately, in a-quartz, the GW 
corrections have been computed only for empty states that are close to the conduction band edge 
[M, B5|, thus we can not verify the accuracy of a scissor-operator approximation far from the edge. 



Recently Mo and Ching |l7| have calculated isotropic Si and O A"-XANES and the corresponding 
LDOS using a supercell approach within an ab initio orthogonalized linear combination of atomic 
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orbitals (OLCAO) method. A good agreement has been obtained between experiments and XANES 
calculations. With regards to LDOS calculations, the situation is quite different. They have shown 
that the Si if-edge can be reasonably represented by Si p LDOS, while O p LDOS surprisingly fails to 
reproduce any feature of the experimental spectrum. This result would mean that the radial dipole 
matrix element are strongly energy dependent in that case. This last consideration may be taken 
cautiously since the energy dependence of the dipole matrix element as well as the LDOS are widely 
dependent on the method used for their calculation. In particular, the LDOS depend on the size of 
the integration region used for the angular momentum projection. Only if this region is sufficiently 
small, the LDOS are independent of the size of the region, appart from a multiplicative factor. We 
calculated Si and O p LDOS, using an integration sphere of 0.4 A radius. The resulting LDOS are 
similar to the corresponding total cross sections computed with our approach, in contradiction with 
the conclusions of Mo and Ching [jf7|. To explain this discrepancy, it should be noticed that Mo 
and Ching used the Mulliken population to compute the LDOS, which corresponds to a much larger 
effective integration region. 

4 conclusion 

We have presented an ab initio framework based on pseudopotentials for calculating XANES spectra. 
The scheme uses self-consistent full-potentials. It was applied to the calculation of K-edge spectra in 
diamond and a-quartz. The results were compared with isotropic (diamond) and dichroic (a-quartz) 
experimental data and a good agreement was obtained. In the case of diamond, our 250 atom su- 
percell calculated spectrum is comparable with the BS calculated one of Soininen and Shirley. These 
first applications have shown that the introduction of the recursion method permits the treatment 
of large supercells. This opens new fields of applications like surfaces, small agregates or amorphous 
materials, for which large supercells are required. The continued fraction may also permit the investi- 
gation of more-extended x-ray absorption regions, without adding significant computation time. The 
results suggest that our scheme is quite robust and could be further successfully applied in XANES 
and ELNES spectroscopies. 

We would like to thank G. Tourillon and C. Laffon for having provided polarized O if-edge 
experimental spectra. We also thank J.J. Rehr for encouraging and fruitful discussions. The english 
was polished by S. Joyce. We are grateful to Ch. Brouder for his thorough reading of the manuscript. 
This is an IPGP contribution #XXXX. 
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Figure 1: Calculated C if -edge X ray absorption spectra in diamond for different supercell sizes, compared 
with experimental data (from Ref. p3|). A 286.1 eV shift was applied to the experimental spectrum. 
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Figure 2: Experimental (dashed line) and calculated (solid line) Si if -edge polarized X-ray absorption 
spectra in a-quartz: (a) <j\\ corresponding to e || [001]; (b) a± corresponding to e _L [001]; (c) XNLD or 
(7ii -erj_. A 1840.7 eV shift was applied to experimental data. 
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Figure 3: Experimental (dashed line) and calculated (solid line) K-edge polarized X-ray absorption 
spectra in a-quartz: (a) a\\; (b) a±. A 531.9 eV shift was applied to experimental data. 




Figure 4: Calculated a± at the silicon K-edge (upper graph) and at the oxygen K-edge (lower graph) 
for different cells: a 2x2x2 supercell including one Is core-hole (solid line), a simple cell (hexagonal unit 
cell) also including one Is core-hole (dashed line), and the ground state hexagonal unit cell (dotted line). 



11 



